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Using a variational method, we derive an exact upper bound for one-magnon excitation energy in 
ferromagnetic spinor gases, which limits the quantum corrections to the effective mass of a magnon 
to be positive. We also derive an upper bound for one-magnon excitation energy in lattice systems. 

The results hold for both Bose and Fermi systems in d dimensions as long as the interaction is local 
and invariant under spin rotation. 
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I. INTRODUCTION 

Effective mass is one of the important concepts in 
condensed-matter physics. This appears in various sys¬ 
tems, such as particles in a periodic potential jl|, Fermi 
liquid [ 2 }, and polaron problems |3J. Since the effective 
mass reflects not only the effects of static structures of 
systems but also dynamical processes between particles, 
it is a key quantity to understand the many-body physics. 

Here, we focus on the problem of whether the interac¬ 
tions between particles increase or decrease the effective 
mass compared with the mass of a bare particle. Liquid 
3 He [1] and /-electron systems 0] are examples that the 
effective mass is heavier than the bare mass. Conversely, 
the effective mass is lighter than the bare mass in a high- 
density electron gas jJQ and a polaron system jf|. Thus, 
it is not obvious that the effective mass is increased or 
decreased by quantum many-body effects. 

Ultracold atomic gases are suitable systems for explor¬ 
ing many-body physics owing to their high experimental 
controllability Q. In particular, ultracold atomic gases 
with spin-degrees of freedom have intriguing features 
originating from superfluidity and magnetism, which 
have been extensively explored @, Q since the experi¬ 
mental realization of spinor Bose-Einstein condensates 
(BECs) fToUIej . In the recent experiment performed by 
the Berkeley group E3. the magnon (spin-wave) dis¬ 
persion and its effective mass in the spin-1 ferromag¬ 
netic BEC were precisely measured. The technique of 
high-precision measurement opens up the opportunity to 
study the spin-wave excitations in the spinor gases. Ac¬ 
cording to the mean-field prediction for the magnon ex¬ 
citation, its dispersion is the same as that of the free 
particle dill. This indicates that the effective mass of 
the magnon is the same as the atomic mass at the mean- 
field level. However, the observed effective mass of the 
magnon is heavier than that of the atomic mass IQ- 

In this paper, using a simple variational method, we 
derive a restriction on the effective mass of a magnon 
in the ferromagnetic spinor gases: the effective mass is 
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never decreased by the quantum many-body corrections 
as long as the interactions between particles are local and 
invariant under spin rotation. We also derive the upper 
bound for the magnon dispersion relation in a lattice sys¬ 
tem. These results hold for both Bose and Fermi systems 
in d dimensions. 


II. CONTINUOUS SYSTEMS 

We consider spin-/ bosons or fermions in a d- 
dimensional space. The field operator is denoted by 
where m = —/, —/ + 1, • ■ • , / represents a mag¬ 
netic sublevel. The field operators satisfy the boson or 
fermion canonical commutation relations, 

hMr),^L( T ' , )]<r = 4>m(r)ipl n ,(r') - 

= S m ,m'S(r - r’), (1) 

= [V'm( r ):^m'( r/ )]<T = 0, (2) 

where a = +1 for bosons and a = —1 for 

fermions. We define the local-density operator 
h(r) = J2 m the local magnetization op¬ 

erator F M (r) = Em 1 »(Wm»A( r )i( r ). and the lo¬ 
cal nematic-tensor operator N^(r) = + 

fvf^mn^lnir^nir)/ 2, where (/^) mn (p = x,y, z) is 
the (m, n) component of the spin-/ matrix. The local 
raising and lowering operators are defined by F±(r) = 

where {f±)mn = ( fx±ify)mn■ 

The Hamiltonian for the system is given by 

H = HKin + -HlZ + -ffQZ + -Hint, (3) 

HKin = f drJ2^ln( r ) ( _ ^T7 v2 ) ^™{r ), (4) 

#LZ = -pM z , (5) 

Hqz = qN zz , ( 6 ) 

where M is the atomic mass, p and q are the mag¬ 
nitudes of the linear and quadratic Zeeman energies, 
M a = f drF a (r) (a = x,y,z,+,—) is the total mag¬ 
netization operator, and N a g = J drN a p(r) is the total 
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nematic-tensor operator. The interaction Hamiltonian is 
defined by [|| 



= ^ E C ny^L( r )^L'( r )^n'(r)i) n {r), 

m,n,m' ,n' 

(7) 

Ajr M {r)= E f,nv,f,rn')ip m (r)4>m'(r), (8) 

m,m' 

C%?' = J2 9r </, m; /, rri\E, M) (E, M\f,rv, /, n '), 

T,M 

(9) 

where T = 0, 2, • • • , 2[/J and A4 = —E, —E + 1, • • • , E 
are the total spin and its z component of two collid¬ 
ing particles, respectively, [_‘J is the floor function, qt 
is the coupling constant for the spin-J 7 channel [20(, 
Ajrj^(r) is the pair annihilation operator with the spin 
state (E, M) at position r, and (E, M |/, m; /, in') is the 
Clebsch-Gordan coefficient. The total particle number, 
the total momentum, and the z component of the total 
magnetization are the conserved quantities: 

W,N]=0, [H,P]= 0, [H,M Z ] = 0, (10) 

where N = J drn(r ) and P = 
(—ih/2) f dr S m [^m( r )W’ m (r) - H.c.] (H.c. denotes 
the Hermitian conjugate). 

We consider two Hilbert subspaces dip and %g W . 
Hf is spanned by the eigenvectors of the Hamiltonian 
{|} that satisfy 7V|Tp) = fVlTp), P|Tp) = 0, 
and M z \^f) = fN |Tf). In this state, all N parti¬ 
cles occupy the magnetic sublevel m = /, i.e., the fully 
spin-polarized ferromagnetic state. Hg W is spanned by 
the eigenvectors of the Hamiltonian {I’I'gvv)} that sat¬ 
isfy N\^ w ) = N |^| w ), P|^sw) = fifc l^sw). and 
Mjl^gw) = (fN — 1)|tHg W ), where k is the wave num¬ 
ber of a spin wave. In "Hg W , N — 1 particles occupy 
the magnetic sublevel m = /, and one particle occupies 
m = f — 1. 

We assume that the ground state |g.s.) is the fully fer¬ 
romagnetic state, i.e., |g.s.) is the lowest-energy eigen¬ 
state in Hf, where |g.s.) is normalized as (g.s.jg.s.) = 1. 
The ground-state energy is denoted by E 0 . We de¬ 
fine the one-magnon excited state |e.s.) as the lowest- 
energy eigenstate in 'Hg W , where |e.s.) is normalized as 
(e.s.je.s.) = 1. The eigenenergy of the one-magnon ex¬ 
cited state is denoted by e| w + Eq, where e| w >0. A 
similar definition of the spin-wave excitation is used in 
Refs. 0,123. 

A. Magnon excitation at the nean-fleld level 

Before presenting our results, we show that the magnon 
excitation energy for spin-/ bosons is the same as the dis¬ 


persion of a free particle in the mean-field approximation. 
For simplicity, we consider the case of p = q = 0. The 
Heisenberg equation for the field operator is given by 

Q „ /j 2 

ih—il) m (r,t ) = -^V 2 V> m (r,t) 

+ E C n^’^L'( r ^n'(r,t)^ n (r,t). (11) 

m' ,71,71/ 

In the mean-field approximation jl 8 l , [13, the system 
is described by the Gross-Pitaevskii equation, which is 
determined by replacing the field operator ip m (r,t) in 
Eq. (fill) with the condensate wave function 4' m (r,f): 

d H 2 

ih—^ m (r,t) = -—V 2 ^ m (r,t) 

+ E CZr'** m ,(r,t)* n ,(r,t)* n (r,t). (12) 

m' ,n,n' 

The ferromagnetic solution is given by T m (r,t) = 
e~ l ^ t / h Smj\/no, where /r = Cjjno = 52 /uo is the chem¬ 
ical potential and no is the particle density. Substituting 
e _I/it / s [y / no<5 m j + dT m (r,t)] into Eq. (fT21) and neglect¬ 
ing the higher-order terms of <5'P m (r,f), we obtain the 
equation of motion for the m = / — 1 component: 

^ t ^f-i(r,t) = -^V 2 <5^/_i(r,t), (13) 

where we used Cjj 1 ^ = 0 and = 

< 72 / / 2 - This equation yields the magnon excitation en- 
ergy e^ IF = h 2 k 2 /2M , which is the same as the dispersion 
relation of a free particle. 


B. Result 1 (continuous systems) 

The one-magnon excitation energy e| w satisfies the 
following inequality: 

P ~ (2/ - 1 )q <el W <e° k +p- (2/ - l)g, (14) 

where f° k = H 2 k 2 /2M is the free-particle energy. 

In the long-wavelength limit, we can expand e| w as 
e fc W = P — (2/ — l)g + H 2 k 2 /2M* + 0(|fc| 3 ), where M* 
is the effective mass of the magnon. Therefore, Eq. m 
gives the inequality 

M* > M. (15) 

Thus, the effective mass never becomes lighter than the 
bare atomic mass. Since the upper bound in Eq. m 
coincides with the excitation energy of the noninteract¬ 
ing systems, the contact interaction always enhances the 
effective mass of the magnon. 
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C. Proof of result 1 


III. LATTICE SYSTEMS 


Our strategy for proving Eq. m is based on a varia¬ 
tional method. We note that a similar method was used 
in Ref. [H] for two-component Bose systems with SU(2) 
symmetry. For |'I') £ "Hg W , we can show that the follow¬ 
ing inequality holds: 


<*!*> 


> e SW 
— £ k 


+ Eq. 


(16) 


Here, we take |H/) = E_(—fe) |g.s.) as a trial wave func¬ 
tion, where E_(fc) = J dre~ lkr F_(r) = F+(—k). One 
can easily show that F_(— k) |g.s.) £ 'Hgw- Using the re¬ 
lations H |g.s.) = Eq |g.s.) and F+(k) |g.s.) = 0, Eq. (fl6l) 
is rewritten as 

C SW < (g-s-l [U+( fc )> [H,F-(-k)]] |g.s.) 

* - (g.s.|[F + (fc),F_(-fc)]|g.s.) 

To calculate the numerator of Eq. 113, we first con¬ 
sider the case of p = q = 0. In this case, the equation of 
continuity for the magnetization is given by [241 ] 

j^,(r,i) + V- Jf n (r,t) = 0 (18) 

in the Heisenberg representation, where J® pm (r) = 

-ih/(2M) E rai J(W-AWVi(»')-H.c.] is the spin- 
current operator. As in the case of the derivation of the 
/-sum rule Q, the equation of continuity in the k space 
and the Heisenberg equation for F_(k,t) gives 

[H,F_(-k)} = Hk ■ Jl pin (-fc), (19) 


where Jl pin (fc) = / dre~ ikr [Jf m {r) - iJ^ in (r)\. The 
detailed calculations of the interaction terms are shown 
in the Appendix [A] After some algebra, we obtain the 
numerator of Eq. as 

(g.s.| [F + (k),[H,F_(-k)]\ |g.s.) = 2 fNel (20) 


where we used M f dr Jf pin (r) |g.s.) = fP |g.s.) = 0 and 
N zz |g.s.) = f 2 N\g.s.). The denominator of Eq. 1171) 
becomes 


<g.s.|[F + (fe),F_(-fe)]|g.s.) = 2/IV. (21) 


We thus obtain the upper and lower bounds for p = q = 
0: 0 < e® w < e°. 

The case of p ^ 0 and q ^ 0 can be derived as fol¬ 
lows. When Hlz and Hqz act on the state vector be¬ 
longing to the Hilbert subspace Hf or "Hg W , they re¬ 
duce to —pfN and qf 2 N for Tip and —p(fN — 1) and 
q(f 2 N — 2/ + 1) for Rg W - Therefore, the linear and 
the quadratic Zeeman terms lift only the eigenenergy. 
The increase in the eigenenergy is given by (e.s.| Hlz + 
Hqz |e.s.) — (g.s.| Hlz + Hqz |g.s.) = p-(2f-l)q. There¬ 
fore, we obtain the inequality m- The result in the pres¬ 
ence of the uniform magnetic field is due to the fact that 
a magnon in the magnetic field is the massive Nambu- 
Goldstone mode [25l - t27| . 


We next consider the spin-/ bosons or fermions in 
d-dimensional lattice systems. The annihilation (cre¬ 
ation) operator at the lattice site R is denoted by 
aR.m («H m )> which satisfies the commutation relations 
[n.R,mi d R , m /]& and [dR^ mi q.r' 

= 0. We define the local density op¬ 
erator h(R) = the local magnetization 

operator F a (R) = Em,n(/a)™®km“ i^ ’" , and the lo ~ 
cal nematic-tensor operator N a p(R) = + 

fpfa)mn^ R m a,R^ n /2. We also introduce the annihilation 

(creation) operator of the Bloch state bk <m (bl m ), where 
hk is the crystal momentum. The relation between d Rtm 

and b k ,m is given by d R , m = (l/y/Np) Y,k e * fc R ^k,m, 
where N-p is the number of sites. 

The single-band Hubbard Hamiltonian of the system 
is given by 

H = Hkiu + Hfz + Hqz + Hint, (22) 

Hk in = tR t R'a Rrn a,R' e {k)b krn bk,m, 

R,R r ,m k,m 

(23) 

Hlz = —pM z , (24) 

Hqz = qN zz , (25) 

where t RR > = t R > R £ K is the hopping parameter de¬ 
pending only on the distance \R — R'\, e(fc) is the bare 
band dispersion, M a = Yr F a (R) is the total magne¬ 
tization operator, and N a p = YjR^ocp{R) is the total 
nematic-tensor operator. The interaction Hamiltonian is 
defined by 

Hint — £ ^A^ m (R)Ar M (R), (26) 

R,M,R 

where gjr is the interaction strength of the spin-J 7 channel 
and Ar M (R) = (T 7 , M\ f, m; /, m') a R ,mQ'R,m' IS 

the pair annihilation operator of the spin state (J 7 , A4) 
at site R. We note that the total particle number op¬ 
erator N = Y^r^R)’ the z component of the total 
magnetization operator M z , and the translation opera¬ 
tor Tr commute with H , where Tr satisfies the relation 

TrCIR' , m T R = d R + R ' ,m. 

In the same manner as in the continuous systems, we 
use the two Hilbert subspaces Hf and %g W . Hp and 
■Hg W are spanned by the eigenvectors of the Hamilto¬ 
nian I'I'f) and |4>g W ) that satisfy 7V|Tf) = A^I^f), 
M z |*f> = /^|T F ),andT jR |^ F ) = |^f) and iV|^ w ) = 

N\*y, Msl’I'sw) = (fH— l)|vE'g W ), and T R \^ k w ) = 
e~ lklt \^ k w ), respectively. We assume that the ground 
state |g.s.) is the lowest-energy eigenstate in Hp, nor¬ 
malized as (g.s.|g.s.) = 1, and its energy is Eq. We also 
assume that the one-magnon excited state is the lowest- 
energy eigenstate in H k w with an energy e| w + Eq. 
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A. Result 2 (lattice systems) 


The one-magnon excitation energy satisfies 
V ~ (2/ - 1 )q < e| W < e(fe) + p - (2/ - l)g, (27) 
where 

e( fc ) = + fe ) - £ ( fc ')] (g-s-| rife-,/ |g.s.), (28) 

k' 


with fik.m = b k In the long-wavelength limit, e(k) 

can be expanded as 

e-w*s < 29 ) 


„ 2W v 


[m--(fc)] * = 


i a 2 

ft 2 dkidkj 


~( k ), 


(30) 


where we assume the space-inversion symmetry of the 
system and ) is an effective mass tensor. 

The upper bound in Eq. m reflects the bare band dis¬ 
persion, in contrast to the uniform continuous systems. 
This is due to the lack of the continuous translational 
symmetry of the system. We note that a similar situa¬ 
tion appears in the /-sum rule for lattice systems [2§j. 


B. Proof of result 2 

The proof is almost the same as that for the continuous 
systems. In Eq. urn we take IT) = F_(—fc) |g.s.}, where 
F-(k) = ^2 R e~ tk ' R F-(R). The upper bound of the 
lattice systems is given by the same expression as the 

continuous systems urn 

As in the case of the continuous systems, we first con¬ 
sider the zero magnetic field p = q = 0. The commutator 
[H,F-(—k)] and its double commutator are calculated 
to be 

[H,F_(-k)j 

= (f-)mn[e(k' + k)-e(k')}b{, +km b k ^ n , (31) 

k' ,m,n 

[F + (k),[H,F_(-k)]\ 

= /(/ + !) 'y [ — k) + e(k' + k) — 2e(k')]n k ',m 

k' ,m 

- ifz)mn[e(k' -k) + e{k' + k) - 2e{k'))b{, m b k ^ n 

k' ,m,n 

+ J2 (f*)mn[e(k' + fc) - e(fc' - k)]bl m b k , tn . (32) 

k' ,m,n 

Using the expectation values for the ground 
state, (g-s.| n fc , m |g.s.) = S mJ (GS| h kJ |g.s.) 

and (g.s.| (, fl)mnb{'J) k ,n |g.s.) 

f s SmjS n j (g-S.| bl f b k j |g.s.) (s = 1, 2), we obtain 

(g.s.| [F+(k),[H,F-(-k)]} |g.s.) = 2 fN~e(k). (33) 


The denominator in Eq. C3 is 
(GS| [F + (k),F-{-k)] |GS) = 2/A, which is the 

same as that of the continuous systems. 

The upper and lower bounds for lattice systems at the 
zero magnetic field thus become 0 < e| w < e(fc). In a 
manner similar to that in the continuous systems, the 
magnetic field lifts only the eigenenergy. Therefore, we 
obtain inequality (l27l) . 


IV. CONCLUDING REMARKS 

We have derived the exact upper bound of the one- 
magnon excitation energy for the ferromagnetic spinor 
Bose and Fermi gases in both continuous and lattice sys¬ 
tems, respectively. This result holds for arbitrary spin 
and spatial dimensions under the following assumptions: 
(i) the ground state is fully spin polarized and has zero 
(crystal) momentum, (ii) the Hamiltonian has continu¬ 
ous (discrete) translational symmetry in continuous (lat¬ 
tice) systems, and (iii) two-body interaction is the s-wave 
contact one and is isotropic in spin-space. Since our re¬ 
sults are based on the local spin-conservation law, anal¬ 
ogous to the /-sum rule due to the particle number con¬ 
servation |2:J, the inequality (fT4l) is one of the universal 
properties of the ferromagnetic spinor gases. Inequalities 
(m and m serve as a check for the validity of quan¬ 
tum many-body corrections. We note that the quantum 
many-body correction to the magnon mass recently cal¬ 
culated in Ref. (2(> satisfies inequality (HU). 

From inequality (HU), we can conclude that the effec¬ 
tive mass of a magnon never becomes lighter than the 
atomic mass. Since the equality M* = M holds for non¬ 
interacting systems, we find that the contact interactions 
always increase the effective mass of a magnon. If one ob¬ 
serves the effective mass of a magnon that is lighter than 
the atomic mass, its origin is not the isotropic contact 
interaction but other effects, such as the magnetic dipole- 
dipole interaction, which is long ranged and anisotropic. 
In fact, we showed that the effective mass of a magnon 
can decrease due to the magnetic dipole-dipole interac¬ 
tion [30 . 

Our results are valid not only for ferromagnetic spinor 
Bose gases but also for ferromagnetic Fermi gases. In 
contrast to the spinor Bose gases (see Ref. [3lj). the ferro¬ 
magnetism in Fermi systems is a generally difficult prob¬ 
lem since the ferromagnetic phase is realized in a strongly 
correlated regime. The Nagaoka state is an example of a 
fully spin-polarized ground state in the spin-1/2 Hubbard 
model [iij (33 and SU(A) Hubbard model [34]. For two- 
component (/ = 1/2) Fermi gases, the Stoner instability 
does not occur due to the competition with the pairing 
instability [35j . Fermi gases with higher spin have already 
been realized in the experiments [3(1 [37]. The ferromag¬ 
netic phase transition in the SU(A) Fermi systems has 
been theoretically discussed [3^ ]. 

Future work is to improve the upper and lower bounds 
of the excitation energy. Our upper bound (ED does 
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not depend on the interaction strength. An appropriate 
modification of the trial wave function allows us to incor¬ 
porate the effect of interactions into the upper and lower 
bounds, which may corroborate the Beliaev theory {29| . 
Another prospect will be the upper and lower bounds for 
the spin-wave excitations in other magnetic phases such 
as polar and antiferromagnetic phases. 
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Appendix A: Calculation for Commutators 


I 

In this appendix, we present the detailed calculations for the commutation relation [Hi nt , F^(r)] (/1 = x,y,z). 
According to Ref. 39j, the interaction coefficient can be written as 


/Oh 

n 


= An 


L/J 

E 

k =1 vi, 


,n' H"~ 'y ' y ' ^k{fvi'' m fu k )mn{fvi' m 'fv k )r 


(Al) 


where Aq and A& are the linear combination of gj The interaction Hamiltonian (0 reduces to 


Hint = y fd,r: n{r) I 2 : + ^ E f dr ^ ,u„ (t 

^ fc= 1 d Vl,— ,V k 


L/J 


_ A(0) V-' M k > 

= ^Int T n rrlt ’ 


fc= 1 


MiL, k (r) = ! • • • 


(A2) 

(A3) 

(A4) 


where 

HRnt > 


: : denotes the normal ordering and NvJ.. Vk 

Kir)]: 


(r) is the rank-fc local nematic-tensor operator. First, we calculate 


A, 


[HinlK( r )] = -y E / dryf^rn'nKLiryKnir'Knir’KmiryKln'irKn'ir)] 


771,71,771' ,7l' ' 


An 


= 0. 


E (U)m’n’ -adn,n4L'{ r )^lni r )Kir)yprn{r) - Sm,n’4>ln' ( r )^n( r )'Mr)^ m (r) 


771,71,771’ ,71' 


+6 rn , m KL(r)'<p} l (r)'>Pn(r)’4>n'(r) + a5 n , m KLi r )Kn(r)4>m(r)4>n'('> 


(A5) 


Next, we calculate [i/^, A /t (r)]. We use a short-hand notation, 


( Tk)mn — ifu i ‘ ‘ ‘ fi/ k )r 


(A6) 
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The commutator [H^, F^r)] is calculated to be 

= j dr' £ £ (T k ) mn (T k ) m , n ,(f^ m „ nll 

m,n,m',n',m",n" v\,--- ,Vk 

X [4>L{ r ')iL' ( r ') Vv {r')^n (r% (r)i> n » (r)] 


= ^2 ^2 (Tk)mn{Tk)m'n' (fu)m"n" 

m,n,m' ,n' ,m" ,n" z>i, • - - , 

x -<r8 m ' in >'$2>( r )$L( r )$n>(r)i>n(r) ~ 5m,n"^L'( r )^m'( r )^n'( r )V i 2 3 4 5 6 ™(r') 
+5„, m '4L( r )$L'( r )$n'{r)$n"{r) +a6 n > !m ^} n {r)^2( r )'4’n(r)'ipn''(r 
= Afe H H (rfc)mn[Tfc,/, i ]m'n''0m(»’)V i L( T ')^™'( r )^( r )- 

m,n,m' ,n' z-'i,-" i^fc 


(A7) 


Here, the commutator [T/, f t ,] rnn becomes 


[Tk, ffi]mn — [fv 1 fv 2 ‘ ‘ ' /i/fc > //j]r 


— ‘ ' ‘ ^i-1 [/^i f v k)rrm 

j =1 
k 

= 'y 1 ^ 1 (/in ■ ■ ■ fvj-ifxfvj+i ■ ■ ■ fvk ) 7 


(A8) 


1=1 A 

where e^A is the completely antisymmetric tensor. From Eqs. (IA7I) and (I A8 [) , we obtain 

k 

[-®Int > -^A‘( t ’)] = Afc • • • /i/j • • • /i/Jmij(Ai ‘ ‘ ‘ /a ‘ ‘ ‘ fv k )m'’ 

m,n,m',n' ,i/fc,A,7=l 

X ^(r)^,(r)^„/(r)^ n (r) 
k 

= ^ ^ ( ^*Ei/j/ia(/ki"'/i/j " '/i/Jmn(/i/i '"/A'"/i/Jr 


171,71,171' ,7l' V 1 ,••• ,ZZfc ,A J = 1 




= 0. 


(A9) 


Therefore, the commutation relation [i/i nt ,-FL(r)] = 0 holds for spin-/ bosons and fermions with the contact 
interaction. In the same manner, we can show [Hint, F/(i?)] = 0 for the lattice system. 
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